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SUMMARY 

A two-dimensional computational code, RPLUS2D, which was developed for the reactive 
propulsive flows of ramjets and scramjets, was validated for two-dimensional shock-wave/turbulent- 
boundary-layer interactions. The problem of compression comers at supersonic speeds was solved using 
the RPLUS2D code. To validate the RPLUS2D code for hypersonic speeds, it was applied to a realistic 
hypersonic inlet geometry. Both the Baldwin-Lomax and the Chien two-equation turbulence models were 
used. Computational results showed that the RPLUS2D code compared very well with experimentally 
obtained data for supersonic compression comer flows, except in the case of large separated flows 
resulting from the interactions between the shock wave and turbulent boundary layer. The computational 
results compared well with the experiment results in a hypersonic NASA P8 inlet case, with the Chien 
two-equation turbulence model performing better than the Baldwin-Lomax model. 


INTRODUCTION 

The performance of high-speed inlets is often significantly affected by the interaction between the 
shock wave and turbulent boundary layer. Shock-wave/boundary-layer interaction is one of the most 
difficult problems to solve numerically. The purpose of the present investigation is to test the capability of 
the RPLUS2D code to predict supersonic compression comers, and then to apply the code to a realistic 
high-speed inlet geometry. 

The problem of supersonic compression comers was studied experimentally by Settles et al. (refs. 1 
and 2), who obtained experiment data for two-dimensional compression comers of 8°, 16°, 20°, and 24°. 
The incoming boundary layer had an edge Mach number of 2.85 and a Reynolds number (Re g ) of 
1.7xl0 6 /m based on boundary layer thickness. 

The numerical simulation of supersonic compression comers has been carried out by many 
researchers (refs. 3 to 7), but success was limited to attached flow and with small separated flows. The 
accurate prediction of large, shock-induced separated flows is still a challenge. Although turbulence 
modeling in shock-wave/boundary-layer interaction has been the subject of many investigations (refs. 6 
and 7) during the past few years, there is still considerable doubt as to the general validity of any 
particular turbulence model. Keeping this in mind, the RPLUS2D code was applied in various situations 
that ranged from a nonseparated (i.e., attached) flow situation (8° ramp) to a large, shock-induced 
separated flow (24° ramp). 

The RPLUS2D code solves the full two-dimensional, Reynolds-averaged, Navier-Stokes equations. 
The Baldwin-Lomax (ref. 8) and Chien two-equation (ref. 9) turbulence models were used in this study. 
This report compares the computations performed for M q = 2.85, where a = 8°, 16°, 20°, and 24°, and 
where Re Q = 7.3xl0 7 /m with the experiment data of Settles et al. (ref. 1). 


♦National Research Council — NASA Research Associate at NASA Lewis Research Center. 
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Cruise conditions of the NASA P8 inlet (ref. 10) are typical of a hypersonic air-breathing vehicle 
and so it was selected as a test case of hypersonic inlet flows for this study. The free-stream Mach number 
was 7.4, the free-stream unit Reynolds number (Re o ) was 8.86xl0 6 /m, and the free-stream total 
■temperature (Tt ) was 8 1 1 K. Wedge and cowl surfaces were cooled to provide a uniform surface 
temperature of 0.375 Tt Q . The NASA P8 inlet has been studied numerically by many investigators 
(refs. 1 1 to 13). Recently, Kapoor et al. (ref. 14) compared the abilities of various turbulence models to 
predict hypersonic inlet flows when applied to the NASA P8 inlet. The geometry of the inlet model tested 
is shown in figure 1 . The inlet model was a Mach 7.4, rectangular, mixed-compression design with exiting 
supersonic flow and had an internal compression ratio of 8; it is referred to as the P8 inlet. The purpose of 
the present investigation was to test the capability of the RPLUS2D code to predict hypersonic inlet 
flows. 


SYMBOLS 

M Mach number 

p static pressure 

P p Pitot pressure 

Pt total pressure 

Re Reynolds number 

Tt total temperature 

u t friction velocity, 

X axial distance from the leading edge of the centerbody 

XREF inlet cowl height 
y + law-of-the-wall coordinate, u Y/v 

J T (0 

Y vertical distance from the centerbody 

a ramp angle 

8 boundary layer thickness 

v kinematic viscosity 

p density 

T shear stress 
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Subscripts: 


o tunnel free-stream condition 

T shear stress 

to evaluated at wall 
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RPLUS2D CODE 

The RPLUS2D code was developed to study mixing and chemical reaction in the flowfield of 
ramjets and scramjets (refs. 15 and 16). It employs both an implicit finite volume and a lower, upper- 
symmetric, successive overrelaxation (LU-SSOR) scheme which solves the full two-dimensional, Reynolds- 
averaged, Navier-Stokes equations and the species transport equations in a fully coupled manner. Yoon 
and Jameson (refs. 17 and 18) initially developed the LU-SSOR scheme for nonreacting flows and 
used extensive testing to show that it was very robust and efficient for transonic and supersonic flows. 

We feel that the RPLUS2D code has the potential to provide a substantial speed advantage over the 
Navier-Stokes codes in current use. 

A switching parameter in the RPLUS2D code allows it to be used for either a reacting or nonreacting 
flow. In this study, air was treated as a single species, nonreacting gas. 


Supersonic Compression Comer Flows 

Experimental background, the method used for computational solution, computational results, and a 
discussion of the results of this test of the RPLUS2D code capability to predict supersonic compression 
comers are included herein. 

Experimental background . — Experiments were conducted in the Princeton University 20 by 20 cm- 
high Reynolds-number, supersonic wind tunnel. Compression comer models of 8°, 16°, 20°, and 24° were 
tested on the wall of the wind tunnel (ref. 1). The uniform free-stream conditions were: Mach number of 
2.85, stagnation pressure of 6.8 atm, stagnation temperature of 268 K, and average free-stream unit 
Reynolds number of 7.3x1 0 7 /m. 

The incoming turbulent boundary layer for the four compression-comer experiments had an overall 
thickness (8) of about 2.3 cm. 

Method of solution . — The computational grid used in this study was 151 by 91. The grid was 
uniform in the X direction. To resolve the viscous layer, the grid lines were packed close to the ramp wall 
with hyperbolic tangent functions such that the first grid line was located at ay + of approximately 5.0 
away from the wall. 

The incoming boundary-layer profile was approximated by calculating the development of a 
turbulent boundary layer on a flat plate to the point where the boundary-layer thickness was equal to the 
experimental value. 

The computations were performed on the CRAY-YMP supercomputer at NASA Lewis Research 
Center. The RPLUS2D code typically required a total of 10 min of CPU time to achieve global 
convergence for the Baldwin-Lomax model, which increased by approximately 20 percent when the 
Chien two-equation model was used. 

Results and discussion . — Computational results are presented in the form of surface pressure 
distributions on the ramps. Skin friction calculations carried out previously by Lee (ref. 1 9) were shown 
to be consistent with the experiment data. Therefore, only ramp surface pressure distributions are 
presented in this report. 
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Figure 2(a) shows the surface pressure distribution for an 8° ramp for both the Baldwin-Lomax and 
the Chien two-equation turbulence models. As seen in the figure, the numerical results compare very well 
with the experiment data. Both turbulence models predict almost identical results. 

Figure 2(b) presents the surface pressure distribution for a 16° ramp. The experimentally obtained 
results show a small region of separated flow just upstream of the comer. The computational results 
compare fairly well with the experiment data. Both turbulence models predicted nearly similar results. 

Surface pressure distribution for a 20° ramp angle is illustrated in figure 2(c). The experimentally 
derived results show significant flow separation with the length of the separated region extended for about 
65 percent of the incoming boundary layer thickness. The computational results compare fairly well with 
the experiment data. The Chien turbulence model performed better than the Baldwin-Lomax model. 

Surface pressure distribution for a 24° ramp compression comer is outlined in figure 2(d). The 
results of the experiment show a massive flow separation with the extent of the separated region slightly 
more than twice the incoming boundary layer thickness. The computational results did not compare very 
well with the experiment data. Both the Baldwin-Lomax and Chien two-equation turbulence models 
failed to predict shock-induced, large separated flows. 

It is evident from this study that the RPLUS2D code results compare reasonably well with the 
experiment data, except in the case of large separated flows resulting from the shock-wave/boundary- 
layer interaction. It is still a challenge to accurately predict such flows; more effort is required to develop 
better turbulence models suitable for this problem. 

Hypersonic Inlet Flows 

Included in this section of the paper are the experimental background, method used for 
computational solution, computational results, and a discussion of the results of the RPLUS2D code 
application to a realistic hypersonic inlet flow. 

Experimental background . — An experimental investigation was conducted at the NASA Ames 3.5-ft 
Hypersonic Wind Tunnel to determine the internal flow characteristics of a typical inlet on a hypersonic 
air-breathing vehicle operating at cruise conditions. The P8 inlet cowl height (XREF) was 18.33 cm, fore- 
body length was 81.28 cm, overall length was 136.2 cm, and width was 35.56 cm. A 6.5° forebody wedge 
was designed to produce an oblique shock wave that passed just outside the cowl lip and delivered flow at 
Mach 6.0 at the entrance of the inlet. The wedge was cooled to provide a uniform surface temperature of 
0.357 Tt , where the free-stream total temperature (Tt Q ) was 81 1 K. The free-stream unit Reynolds 
number (Re o ) was 8.86x1 0 6 /m. The experimentally determined laminar-turbulent transition point on the 
centerbody was approximately 40 percent of the distance between the wedge leading edge and the inlet 
entrance. The laminar-turbulent transition point in the cowl boundary layer was approximately halfway 
between the cowl leading edge and the throat station. Details of the experiments can be found in 
reference 10. 

This realistic inlet geometry, with strong viscous-inviscid interactions and the availability of 
extensive experiment data, provided an excellent opportunity to verify the ability of the numerical 
algorithm and turbulence models to predict a hypersonic inlet flow field. 

Method of solution . — The computational grid used in this study was 221 by 91 and nonuniform in 
the X direction. The grid was packed on both ends from the wedge leading edge to the cowl leading edge 
and was also geometrically stretched from the cowl leading edge to the outflow boundary. To resolve the 
viscous layer, the grid lines were packed close to the centerbody and cowl walls using hyperbolic tangent 
functions such that the first grid line was located at a y + of approximately 4 0 away from the walls. 

The laminar-turbulent transition points on the ramp and cowl surfaces for the Baldwin-Lomax model 
were manually simulated from experiments by setting eddy viscosity equal to zero for laminar flows. 

The computations were performed on the CRAY-YMP supercomputer at NASA Lewis Research 
Center. The RPLUS2D code typically required a total of 15 min of CPU time to achieve global 
convergence for the Baldwin-Lomax model, which increased by approximately 20 percent when the 
Chien two-equation model was used. 
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Results and discussion . — The performance of a hypersonic inlet is significantly affected by the 
interaction of the shock waves and turbulent boundary layers. The interaction between the inlet cowl-lip 
shock and centerbody turbulent-boundary-layer, in particular, requires careful analysis since the 
centerbody contours are often designed to cancel the cowl-lip shock. The internal contours of the P8 inlet 
model were designed to provide cancellation of the cowl shock at the centerbody and an isentropic 
compression to the throat. Nevertheless, a reflected shock was found experimentally, which further 
interacted with the cowl and centerbody boundary layers upstream of the throat. This significantly 
affected the flow structure at the inlet throat. The present computations were able to capture these flow 
characteristics. 

The computed density, pressure, and Mach number contours for the Chien turbulence model are 
shown in figure 3. The contours obtained from the Baldwin-Lomax turbulence model were essentially 
similar. Figure 3 clearly illustrates that the cowl shock, after interacting with the centerbody turbulent- 
boundary-layer, was reflected downstream and interacted with the cowl turbulent-boundary-layer. The 
cowl shock was then reflected from the cowl and interacted again with the centerbody turbulent- 
boundary-layer before it left the inlet. 

Further computed results are presented in the form of surface pressure distributions on the 
centerbody and cowl of the inlet model, and the Pitot pressure and total temperature distributions at many 
stations from the inlet entrance to the throat of the inlet. The results are compared with the corresponding 
experiment data from reference 10. 

The surface pressure distributions on the centerbody are shown in figure 4. The axial distances are 
nondimensionalized with the inlet cowl height. The surface pressure distributions indicate that the 
computed results predict the location of the interaction of the cowl shock wave with the centerbody 
turbulent-boundary-layer slightly upstream as compared with the experiment. 

The results obtained using the Baldwin-Lomax model show the existence of a separation bubble on 
the centerbody in the immediate region of cowl-lip-shock/centerbody-boundary-layer interaction. The 
existence of the separation bubble was confirmed by the presence of negative velocities in the separated 
region as shown in the insert of figure 4. No separation was reported in the experiment. Kapoor et al. 

(ref. 14) also reported the presence of a separation bubble when they used the Baldwin-Lomax model in 
the PARC2D code. On the other hand, the Chien model does not show the presence of a separation bubble 
and is able to successfully simulate the complex flow field resulting from the interaction between the 
cowl shock and the centerbody turbulent-boundary-layer. 

The pressure distributions on the cowl surface are presented in figure 5. The computed results show 
the flow compression and impingement of the reflected shock wave on the cowl. The expansion ahead of 
the reflected shock impingement, which is a feature associated with the shock-wave/boundary-layer 
interaction on the centerbody, is also visible in the computed results. The Chien two-equation turbulence 
model performs better than the Baldwin-Lomax model in predicting cowl surface pressure distributions. 

The computational and experimentally derived Pitot pressures where X/XREF equals 5.67, a station 
upstream of the intersection of the cowl shock with the centerbody, are shown in figure 6. The agreement 
of the computed results with the experiment data is generally good. The steep rise in the Pitot pressure 
where Y/XREF equals 0.15 is due to the presence of the cowl shock. The results from the design analysis 
in reference 10 and the boundary-layer thicknesses on the centerbody and cowl surfaces obtained from the 
experiments are also marked in figure 6. The Chien two-equation model compares fairly well with the 
experiment data. 

All of the present computational results overpredicted the Pitot pressures in the central region of the 
inlet. A recent AGARD report (ref. 13) noted that no one has ever matched the experiment data using the 
tunnel conditions stated in reference 10. The various computational results tend to agree with one another, 
but do not match the experiment data. The AGARD report further suspected that the conditions stated in 
reference 10 were different in some way from the conditions that were actually present in the tunnel. It 
should be noted that because of the high sensitivity of hypersonic flows, even small variations in the 
upstream flow field would lead to larger variations downstream. 
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The total temperature profile where X/XREF equals 5.67 is presented in figure 7. The 
experimentally obtained total temperature ratio is less than 1 .0 because it was not corrected for the 
uncertainties in other measurements (ref. 10). The comparison of the computed results with experiment 
data is generally good. The Chien two-equation model compares well with the experiment data, 
particularly near the centerbody. 

The pitot pressure and total temperature distributions where X/XREF equals 6.09, the intersection 
point of cowl shock with centerbody, are presented in figures 8 and 9, respectively. The computational 
results compare qualitatively well with experiment data. 

The pitot pressure and total temperature distributions where X/XREF equals 6.37, a station just 
downstream of the reflection of the cowl shock from the centerbody, are presented in figures 10 and 11, 
respectively. Figure 10 illustrates that the centerbody boundary layer of the experiment has been 
compressed by the reflecting shock wave and is thinner than at the previous station. The reflected shock 
wave emerges from the boundary layer where Y/XREF equals 0.025, as shown by the break in the curve. 
The computations are able to detect the emerging shock wave, but the magnitude and location of the 
emerging shock differs for both turbulence models. Design analysis from reference 10 is also presented. 
The comparison of total temperatures with experiment data is good, as shown in figure 1 1 . 

The pitot pressures and total temperatures where X/XREF equals 6.65, a station upstream of the inlet 
throat, are shown in figures 12 and 13, respectively. The comparison of pitot pressures and total 
temperatures with the experimentally derived data is qualitatively good. 

The inlet performance was obtained in terms of pitot pressure, total temperature, and Mach number 
distributions at the throat. Figure 14 depicts a large variation of the pitot pressure distribution across the 
throat height due to the presence of the reflected shock wave. Results of the design analysis from 
reference 10 are also shown in figure 14. The total temperature distribution at the throat is shown in 
figure 15 and the Mach number profile at the throat is presented in figure 16. The computed results are 
in fair agreement with the experimentally derived data. 


CONCLUDING REMARKS 

The RPLUS2D code was validated for two-dimensional shock wave and turbulent boundary layer 
interactions. Supersonic flow over compression comers and hypersonic flow in the NASA P8 inlet were 
numerically simulated. Both the Baldwin-Lomax and Chien two-equation turbulence models were 
investigated in this study. 

The computational results show that the RPLUS2D code compares very well with experimentally 
derived data for supersonic compression comer flows for the nonseparated and small separated flows, but 
it fails in the case of large shock-induced separated flows. 

The computational results from the RPLUS2D code compare reasonably well with experimentally 
derived results when applied to the hypersonic NASA P8 inlet. Our conclusion from this study is that the 
Chien two-equation turbulence model performs better than the Baldwin-Lomax model for hypersonic 
inlet flows. 
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Static pressure ratio, p/po Static pressure ratio, p/p 0 
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TtQ = 811 °K 
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Axial distance, X/XREF 

Figure 1 .—Diagram of P8 inlet model. 
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Figure 2. — Comparison of computed and experimentally determined pressure distribution for an inlet ramp in supersonic flow where (a) a = 8°, 
(b) a = 1 6°, (c) a = 20°, and (d) a = 24°. 
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Figure 3. — Typical (a) density, (b) pressure, and (c) Mach number contours for P8 inlet. 
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Figure 4. — Surface pressure distribution on the centerbody and 
velocity distribution in separated flow region. 
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Figure 5. — Surface pressure distribution on the cowl. 



Figure 6. — Pitot pressure distribution where X/XREF equals 5.67. 



Figure 8.— Pitot pressure distribution where X/XREF equals 6.09. 



Figure 9. — Total temperature distribution where X/XREF equals 6.09. 
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Figure 10.— Pitot pressure distribution where X/XREF equals 6.37. Figure 13. — Total temperature distribution where X/XREF equals 6.65. 




Figure 11 . — Total temperature distribution where X/XREF equals 6.37. Figure 14. — Pitot pressure distribution at the throat. 




Figure 1 2.— Pitot pressure distribution where X/XREF equals 6.65. Figure 1 5.— Total temperature distribution at the throat. 
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